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In systems where one coordinate undergoes periodic oscillation, the net displacement in any 
other coordinate over a single period is shown to be given by differentiation of the action integral 
associated with the oscillating coordinate. This result is then used to demonstrate that the action 
integral acts as a Hamiltonian for slow coordinates providing time is scaled to the "tick-time" of the 
oscillating coordinate. Numerous examples, including charged particle drifts and relativistic motion, 
are supplied to illustrate the varied application of these results. 
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Hamiltonian dynamics is almost ubiquitous in physics and describes such varied phenomena as celestial mechanics, 
optics, fluid dynamics, quantum mechanics, and charged particle motion in electromagnetic fields. Guiding center 
theory, an approximation of Hamiltonian dynamics for charged particle motion in magnetic fields, describes the motion 
of the particle's cyclotron-orbit averaged position, or guiding center [l[. The guiding center can be thought of as 
a "quasi-particle" subject to new types of forces and manifesting various drifts. We develop a general model, not 
restricted to charged particle motion, of multi-dimensional systems with a periodic variable and find drifts that cannot 
be calculated using guiding center theory which becomes a limited example of the more general model. The model 
shows that the action integral associated with the oscillatory coordinate acts as an effective Hamiltonian for the 
remaining, slow coordinates providing time is measured in clock cycles of the oscillations. We note that Hamiltonian- 
type aspects of action integrals have been previously discussed in specific situations 0] but without developing a 
general demonstration and relying on the detailed equations of motion in their proofs. The model presented here 
generates a hierarchy of "wheels-within-wheels" Hamiltonian systems such that the action integral associated with 
periodic motion at any level in the hierarchy acts as the Hamiltonian for the next slower periodic motion. 

Consider a 2D time- independent Hamiltonian system H(£, P v ) with an ignorable coordinate 77 and where the 
^-motion is periodic, i.e., £(t + At) = £(t), with no limit on the amplitude of £. P v evolves trivially: P v = 0, but the rj 
evolution is in general nontrivial. The period At can be imagined as a clock tick over which 77 undergoes a net change 
A77. We claim that 

A V = -8J/dP v , (1) 

where 

J(H,P V )= I Pt(H,S,P v )dt (2) 

is the action integral [|[ associated with £. Equation (JTJ means that if J(H, P^) is known, the net change of 77 during 
one period of £ can be calculated without having to consider the potentially complicated form of 77. 
To prove Eq. ([1]), first note that there is no contribution from differentiating the integral's bounds, so 

dJ I dP £ 

The differential of H is 



so 



dP L = _dH/dP 1 

dP n dH/dP^ [ ' 

Using Eq. (|5j) and Hamilton's equations in Eq. ([3]) gives 

dJ f dH/dP v 

dP v J dH/dPs * 

If there are other ignorable coordinates in the system, then suitably adjusted versions of Eq. ([6]) apply to each of 
them. Equation ([1]) generalizes the well-known theorem Q that the period of motion is given by a partial derivative 
of J with respect to H, namely 

At = dJ/dH. (7) 

Equation (J7J resembles Eq. (Q} because (t, —H) form a pair of canonical coordinates in extended phase space, so 
Eq. is a special case of the theorem presented here. The drift, or net time evolution, of 77 is clearly 

Ar, dJ/dPr, 

At dJ/dH ' [ ' 



which generalizes the particle drifts associated with guiding center theory. 
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We now relax the requirement that 77 is ignorable and allow the oscillations to evolve adiabatically. We do so by 
coupling the original Hamiltonian, now denoted as -ffi oc , to an external system H cxt that is otherwise isolated. This 
gives a total Hamiltonian 

H^,P i ,r l ,P 11 )=H loc (^P^r l ,P 11 )+H cxt (r l ,P 11 ). (9) 

We presume the system behaves as follows. First, the ^-oscillation is described by the local Hamiltonian H\ oc in 
which i] and P v play the role of slowly varying parameters: d£/dt — dH\ oc / dP^ and dP^/dt = —dH\ oc / dt;. Second, 
the "parametric" coordinates r\ and P v are described by the total Hamiltonian H, so drj/dt = dH/dP n and dP v /dt = 
—dH/drj. The local and external systems exchange energy, but the total energy E = E\ oc (t) + E ext (t) is conserved 
since the entire system is isolated. J is defined as in Eq. ([2]) but is now also a function of 77. As in Ref. Q, we assume 
it is a good approximation to hold the parametric coordinates 77 and P v fixed while evaluating the £ action integral. 
J is an adibatic invariant and is thus conserved. Furthermore, J depends only on H\ oc , i.e. J = J(H\ oc , 77, P, q ) = 
J{H — H ext (r], P v ), 77, P v ), because H\ QC is sufficient to prescribe the £ dynamics. A proof analogous to that of Eq. (fTJ) 
shows 

dJ . „ dJ , „. 

^ = AP - w v = ~ Av - (10) 

Note that J = J(H — H ex %(r), P v ), 77, P v ) depends on 77 and P v both implicitly through H cx t and also explicitly. 
Accordingly, A77 and AP V have two terms: one term comes from the explicit dependence and is the drift of the 
system; the second term comes from the implicit dependence and is the slow change of H cxt - 

Equations (|10p have the makings of a Hamiltonian system with — J serving as the Hamiltonian. They are precisely 
Hamiltonian as follows. We define discretized derivatives drj/dt = Arj/At and dP„Jdt — AP v /At that capture the 
net rates of change of rj and P v . Upon invocation of a rescaled time r normalized by the ^-period: 



Eqs. (fT0| become 



dr = dt/At, (11) 



Thus, — J is the Hamiltonian for the averaged system provided time is measured in units of At. It should be noted 
that r is the angle variable conjugate to J. 

Alternatively, we can obtain a Hamiltonian for the ^-averaged system by solving J = J{Hi oc , rj, P^) for Hi oc = 
Hioc(J,V> Pr/) which upon inserting in Eq. (9) gives 

H = Hi oc (J, 77, P n ) + H ext (r/, P n ). (13) 

Solution of Eq. (13) for J gives J = J(H, 77, P n ). The differential of J using this latter form is 

Equation (14) determines dH/dP v = — {dJ/dP^) / (dJ/dH) etc., so using Eqs. (10) 



dH 


dJ/dP n 


-A77 


drj 


dPr, ~ 


dJ/dH 


At 


lit 


dH 


dJ/dr, 


AP V 


dP v 


drj 


dJ/dH 


At 


dt 



(15) 
(16) 

Thus, H, written as Eq. (|13p . generates the discretized derivatives. The term H\ oc (J, 77, P^) is an adiabatic potential 
the residue of averaging the periodic ^-motion. For systems approximating a harmonic oscillator, J = 2nHi oc /uj(r], P ri ) 1 
so the adiabatic potential is H\ oc — Jlo/2t:, showing that J acts like an electrostatic charge and w^P^) acts like 
an electrostatic potential. The magnitude of this " J-charge" depends on the amplitude of the ^-oscillation. The use 
of — J(H, rj, P v ) as a Hamiltonian with normalized time r and the use of iJ(J, 77, P v ) with regular time are entirely 
equivalent. Practically, though, there are techniques to evaluate J directly @, so using — J as the Hamiltonian spares 
one from inverting J for H , which might not be analytically feasible. 

We now provide examples illustrating applications. Fig. [1] shows an electron moving in the rz-plane and subjected 
to the magnetic field of a current-carrying wire aligned along the z-axis. The z coordinate corresponds to 77 and is 
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FIG. 1. For planar electron motion outside a current-carrying wire, the axial displacement Az can be derived from the radial 
action variable. 



ignorable; the radial motion is periodic and not ignorable because of the magnetic field gradient. The electron displaces 
itself an axial distance Az with every gyration as shown in Fig. [T] Using the characteristic velocity j3 — fj, Ie/2irm 
Q, the Hamiltonian is 

H= i mv2= ^ + (^^hml, (17) 

2 2m 2m 

where R is an arbitrary constant of integration. J can be evaluated exactly Q using the substitution cos 6 = 
(P z — m(3\n(r/R)) jmv and the integral representation of the modified Bessel function /„ (x) = ir^ 1 J* e xcosS cos(n8)d9 @ 
so that 



J = j> p rdr = 2Trmvr gc I 1 I - j , (18) 

where r gc = Rexp(P z /mf3). r gc , plotted as a dashed line in Fig. [TJ is the radial position at which the z-velocit; 
vanishes, as can be checked from P z = mv z + m/3 In r/R. J generalizes the first adiabatic invariant /i = mv\/2B [lj 
and reduces to 27r(m/e)/i when v <C (3, which for this system is the condition for the guiding center approximation to 
hold @- Az and At can be computed using Eqs. (JTJ) and (J7]) and noting that v = ^/2H/m. The exact drift velocity, 
computed without appealing to the guiding center approximation, is 

Az. h{v/p) A 
At Io{v/p) 

Equation (fT9|) holds for orbits of all energies even when the guiding center approximation fails. The v <C (3 limit of 
Eq. (fTS)) reduces to the grad B drift [l(| of the guiding center approximation. 

Next we show how J can be used as a Hamiltonian to give the magnetic mirror force [10]. For a magnetic field 
mainly in the z direction, the cyclotron motion is essentially harmonic oscillation at the gyrofrequency u) = qB /m in 
the perpendicular direction, so we identify 77 with z and H cyLt with P^/2m. Then 

H loc _ H-P z 2 /2m ^ fnvj/2 

J = 2tt =2tt — =2tt— - — , (20) 

to qo/ra qB/ra 



which equals /i except for the factor 2irm/q. If B depends on z then Eqs. (fT5|) and (| 16[) become 

dz dJ/dP z P z /m qB _ P z 
dt dJ/dH qB/rn m m 



(21) 



dP z _ dJ/dz _ qj dB _ dB 

~dT~dJ/dH ~ ~2^i ~dz~ ~ ~^~dz~' ( ' 

establishing the magnetic mirror force without considering the microscopic motion. 

A slightly different approach retrieves the grad B drift. Suppose B =B z (x)z, so A = A y (x)y with B z — dA y /dx. 
We define the x component of the guiding center as the position x gc where v y vanishes: P y — qA y (x gc ). Setting 
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r] = y, Eq. (|T5|) applied to Eq. (|20l) gives 

dy dJ/dP y f dB z 
dt dJ/dH fI \ dx 



dx 



dP v 



(23) 



We then use dx gc /dP y 
thus becomes 



(qB(x gc )) 1 , obtained by differentiating P y = (7Aj,(a; gc ) with respect to P y . Equation (j2"3")l 



cfy _ / d5 2 
dt qB z \ dx 



(24) 



which is the grad B drift. 

A surprising application arises in relativistic mechanics, where it is found that in crossed electric and magnetic fields 
E = Ex and B = Bz with E < Be a particle's z-velocity is modulated by the cyclotron motion, in contrast to the non- 
relativistic situation where v z is constant and independent of the cyclotron motion. The modulation arises from the 

1/2 

periodic addition and subtraction of the E x B drift to the cyclotron velocity, which modulates 7 = (l — u 2 /c 2 ) 
and hence the particle's effective mass; v z then varies because v z = P z /jm and P z is invariant as z is ignorable. 
Using the relativistic canonical momenta P = rrryv + qA with A =Bxy, the electrostatic potential 4> = —Ex, and 
the relativistic Hamiltonian H = c^J (P 
evaluates to 



qAf 



q<j), it is found that the relativistic x-action J = £ P x dx 



2 i2 



J_ _ {BP y + HE/c 
2n ~ 2q(B 2 - E 2 /c 2 ) 3 / 2 



(H 2 /c 2 -P 2 



p2 

y 



m 2 c 2 



2q(B 2 - E 2 /c 2 y/ 2 



Calculating Az and At using Eqs. ([TJ and gives the independent z-drift velocity 

B 2 - E 2 /c 2 



Vd 



Az 
At 



BEPy/c 2 + B 2 H/c 2 



:P, 



(25) 



(26) 



As shown in Fig. (|2"|). this v z drift has been verified by direct numerical integration of the relativistic equation of 
motion di^ymv) / dt = g(E|vxB) which shows that the modulation of v z is typically spiky as 7 1 for a short 
interval during each cyclotron period and then 7 3> 1 for the remaining fraction of the cyclotron period. Clearly, this 
analysis generalizes to force-drifts by replacing E with F/q. 
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FIG. 2. A particle undergoing relativistic E x B motion in the icy-plane with E = 0.95-Bc and initial momentum P x = 0, 
P y = 0.7mc, and P z = 0.3mc. (a) The z-velocity (from numerical integration of the relativistic equation of motion) is non- 
constant, spiking when the particle's icy-motion slows down so that 7 1 
dashed line, calculated using Eq. ((261) . captures the z-drift motion. 



1. (b) Solid line is the numerically integrated z-position; 



Kepler motion provides a non-relativistic and non-electromagnetic example. The radial action is Q 



J r = -2ir |P | + 2irVmmMG/y/2\H\, (27) 
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where is the conserved angular momentum. Equation ([T]) gives A</> = ±27r depending on the sign of P^, immediately 
proving that bounded Kepler orbits are always closed. 

We now present a purely mechanical system which exhibits the equivalent of "magnetic" mirroring. Consider a 
non-relativistic particle in a long groove where the width of the groove varies with position. The Hamiltonian is 



P'z 



Pi 



X 



(28) 



where y is the distance along the groove and x is the distance across the groove. Presuming that the y-position 
changes slowly relative to the oscillations across the groove (i.e., \a\ and |A| are small compared to k) the y-dependcnt 
frequency of x-oscillation is 



(29) 



We identify H\ oc = P 2 /2m + muj(y) 2 x 2 /2, so the x-action is J = 2nH\ oc / 'uj(y), and Eq. (fT3|) becomes 



H 



2m 2 y 2tt 



J. 



(30) 



Equation (fTB| gives an average force — (J / 2tt) du / dy = —yjKa/2irmu> in the y-direction due to the increase in 
rc-oscillation energy where the groove narrows. This is a restoring force and, if sufficiently strong, can overwhelm 
the contribution from A. A negative A corresponds to a potential hill, and if J — the particle falls down the hill. 
However, if J is sufficiently large and a is positive, the particle does not fall down but instead oscillates about the top 
of the hill! This mechanical analog of a magnetic mirror has been verified by direct numerical integration as shown 
in Fig. [3 




m 0.04 - 



-1.0 -0.5 



0.0 

y 



0.5 



1.0 



FIG. 3. A particle in a thin saddle-like groove can undergo oscillatory motion due to narrowing of the groove, x is the direction 
across the groove, y along the groove, and z the vertical. H is given by Eq. (|28(l with m = k = a = 1 and A = —.01, and the 
particle starts at x — y = with v x = .25 and v y — 1. 



For oscillatory y-motion, Eq. (1301) admits an action integral in the y-direction, which we denote by K, that acts as 
a Hamiltonian for the a;-averaged system. This is a two-tier heirachy of action variables, or a wheel within a wheel. 
For the reduced system, J is a conserved quantity, so we develop an analog of Eq. ([1} : 

dK f dP v (H,J,y) f 1 

1 l -dy = -Ar, (32) 



/ dy/di 



where we use Eq. (fTl)) to evaluate dP y /dJ and Eq. (fT2")l to evaluate dJ/dP y . Since r counts re-cycles, —dK/dJ gives 
the number of x-cycles per y-cycle. If this quantity is a rational number, the trajectory is closed. This is of interest 
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when quantizing the system,_as there is sometimes a one-to-one correspondence between periodic classical trajectories 
and quantum energy levels [Tl| • 
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